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Abstract 

We examine the main properties of relic neutralinos in the context 
of the Minimal Supersymmetric extension of the Standard Model. We 
present a whole set of results obtained for the relic abundance and for 
detection rates with a model constrained according to the latest acceler- 
ator data. Predicted detection rates for relic neutralinos are compared to 
the present experimental bounds. A comparison with results obtained 
with more sophisticated theoretical schemes with unification assump- 
tions at Mgut is presented too. We also briefly discuss some properties 
of a Higgsino-like neutralino which has recently been suggested in con- 
nection with a supersymmetric interpretation of a single event at the 
Fermilab Tevatron. 



1 Introduction 



The neutralino is currently being considered as a favourite candidate for Cold 
Dark Matter (CDM). In fact, the neutralino turns out to be the Lightest Super- 
symmetric Particle (LSP) in extended regions of the supersymmetric parameter 
space and then, if protected by R-parity conservation (that we assume to hold), 
it would have decoupled from the primeval plasma and would be present in 
the Universe as a fossil particle. 



The properties of the neutralino vary quite noticeably as one moves in the 
supersymmetric parameter space: from some regions where the neutralino relic 
abundance £l x h 2 (h is the present Hubble expansion rate in units of 100 Km • 
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s _1 • Mpc -1 ) turns out to be large (and detection rates low) to other regions 
where the other extreme situation of a negligible Q x (and of sizable detection 
rates) occurs. It is remarkable that there are domains of the supersymmetric 
parameter space, where the neutralino can provide a relic abundance in the 
preferred density range for CDMEl 



VcDMh 2 = 0.2 ± 0.1. (1) 

However, even in the case when the neutralino does not contribute signif- 
icantly to the average density of the Universe, experimental searches for relic 
neutralinos would be of the utmost interest for a number of cosmological and 
particle-physics aspects. 

In this paper we examine the main properties of relic neutralinos in the 
context ppf the Minimal Supersymmetric extension of the Standard Model 
(MSSMp. Indeed, this scheme provides a most convenient phenomenologi- 
cal framework at the electroweak scale (Mz), without assuming too strong, 
arbitrary theoretical hypotheses. After a short presentation of the theoreti- 
cal scheme, we report the main results concerning neutralino relic abundance 
and predicted detection rates. A comparison with results obtained with more 
sophisticated theoretical schemes is presented too. 

Furthermore, we also briefly discuss some properties of a Higgsino-like neu- 
tralino which has recently been suggested in connection with a supersymmetric 
interpretation of a single event at the Fcrmilab Tevatron. 



2 The Minimal Supersymmetric Standard Model 

The MSSM is based on the same gauge group as the Standard Model and 
contains the supersymmetric extension of its particle content. The Higgs sector 
is modified as compared to that of the Standard Model, because it requires 
two Higgs doublets Hi and Hi in order to give mass both to down- and up- 
type quarks and to cancel anomalies. After Electro- Weak Symmetry Breaking 
(EWSB), the physical Higgs fields consist of two charged particles and three 
neutral ones: two scalar fields (h and H) and one pseudoscalar (A). The 
Higgs sector is specified at the tree level by two independent parameters: the 
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mass of one of the physical Higgs fields (we will use the mass itia of the A 
boson) and the ratio of the two vacuum expectation values, usually defined as 
tan/3 = v 2 /v 1 =< H 2 > / < Hi >. 

Apart from the Yang-Mills Lagrangian, other characteristic elements of 
the model are the superpotential, which contains all the Yukawa interactions 
and the Higgs-mixing term \iH\H 2l and the soft-breaking Lagrangian, which 
models the breaking of supersymmetry 

-C 8of t = Y, m ^ 2 + {[ A -b h abLaH 1 R b + AiXbQaH 1 D b 

i 

+ A u ab h u ab Q a H 2 U b + h.c] - B[lH\H 2 + h.c.} 

+ ]rM 4 (AA + AA). (2) 

i 

Here the notations are the following: the fa are the scalar fields, the Xi are 
the gaugino fields, Q and L are the doublet squark and slepton fields, respec- 
tively, and U, D and R denote the S[/(2)-singlet fields for the up-squarks, 
down-squarks and sleptons, m, and Mi are the mass parameters of the scalar 
and gaugino fields, respectively. A denotes the trilinear and B the bilinear 
supersymmetry breaking parameters. The Yukawa interactions are described 
by the parameters h, which are related to the masses of the standard fermions 
by the usual expressions, e.g., m t — h t v 2 . 



The tree-level scalar potential for the neutral Higgs fields turns out to be 

V = (Af| ri +/i 2 )|i7 1 | 2 + (M| r2 +^ 2 )|ff 2 | 2 -B^(i7 1J ff 2 +h.c.)+quartic D terms. 

(3) 

When the electro- weak symmetry breaking is triggered by H\ and H 2 acquiring 
vacuum expectation values v\, v 2 with the condition vf + v\ = v 2 = Im^jg 1 , 
the potential Vb takes the standard form of a Higgs potential, provided the 
following contraints are satisfied: 



(4) 



(5) 



Notice that the sign of [i is defined according to the convention of Refs. □. 
We remark that Eqs.(^-^|) are expressed at the tree level. However, in the 
calculations that we will report later, 1-loop corrections to Vq are included a, 

As is manifest from the previous formulae, the theoretical model contains 
an exceedingly large number of parameters. In order to have a viable phe- 
nomenological model we reduce the number of free parameters by introducing 
the following usual assumptions: 

• A33 = = A l 33 = A; all other A &'s are set to zero 

• all squarks are taken as degenerate: = mo, except for 

• all sleptons are taken as degenerate with m ; ~ = mo 

• the gaugino masses are assumed to be renormalised as the gauge cou- 
plings, 

Mi : Mi : M3 = ol\ : ct2 ■ 013; thus in particular Mi = (5/3) tan 2 dw^i- 

We emphasize that all these relations are meant to be taken at the Mz 
scale. 

It is worth noticing that any variation in some of these assumptions, if re- 
quired, can be introduced staightforwardly without implying any other modifi- 
cation in the general scheme. The flexibility of this phenomenological approach 
is not shared by other more theoretically-constrained models. 

In force of the conditions previously discussed, the number of independent 
parameters is reduced to seven. We choose them to be: M2, /J., tan (3, tra, mo, mi, 
A. In terms of these parameters, the supersymmetric space is further con- 
strainted by the following requirements: 



• all experimental bounds on Higgs, neutralino, chargino and sfcrmion 
masses are satisfied (taking into account also the new data from LEP2a) 

• the neutralino is the Lightest Supersymmetric Particle (LSP) (i.e., re- 
gions where gluino or squarks or sleptons are LSP are excluded) 

• constraints due to the b —> s + 7 process are taken into account! 
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• the neutralino relic abundance does not exceed the cosmological bound, 
i.e. Q x h 2 < 1. 

As usual, any neutralino mass-eigenstate is written as a linear superposi- 
tion 



X% = a il + h% + c%H s + diH a (6) 

where y,Z are the photino and zino states and H s ,H a are defined by H s = 
s'm(3Hi + cos/3H%, H a = cos(3H^ — sin/3 H^, in terms of the higgsino fields H^, 
H%, supersymmetric partners of the Higgs fields iff, H%- In the following we 
will simply call neutralino the eigenstate of lowest-mass and use \ = Xi- Figjl] 
gives the isomass contours in the (J, — M2 plane. 

In the same plot, the region excluded by the LEP data and the contour 
lines for the composition parameter P = a\ + b\ are also displayed. 

Various properties of relic neutralinos in the MSSM have been considered 
by a number of authors. Some of the most recent ones are given in Refs&Q. In 
this review, we present a whole set of results concerning the main properties for 
relic neutralinos (relic abundance and detection rates) performed in the same 
consistent theoretical frameworkfl The techniques emckcpd for the calculation 
of these various quantities are those presented in Refs&EjEi The experimental 
constraints have been updated to include the latest LEP2 dataa. Only the 
predictions for the p/p ratio are not given here, since they_are presented and 
discussed in an accompanying paper in these Proceedings 113. Our results for 
the neutralino relic abundance and detection rates are provided in the form of 
scatter plots. The independent parameters of the model have been varied in 
the following ranges: 



10 GeV 


< M 2 


< 500 GeV 


10 GeV 


<M 


< 500 GeV 


50 GeV 


< m A 


< 500 GeV 


100 GeV 


< m 


< 500 GeV 


100 GeV 


< mi 


< 500 GeV 
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-3 < A < +3 

1.1 < tan/3 <55 (7) 

Let us mention here that in more sophisticated supersymmetric models 
the properties at the Mz scale are derived from properties at the Gran Uni- 
fication scale Mqut by employing the renormalization group equations. One 
of the nicest features of this approach is the link between the supersymmetry 
breaking and the EWSB, which is induced radiatively. These supersymmet- 
ric schemes are usually based on assumptions of unification at Mqut, not 
only for the gauge coupling constants, but also for Yukawa couplings and for 
the soft-breaking parameters (scalar masses, gaugino masses, trilinear cou- 
plings). These assumptions have very strong consequences for the neutralino 
phenomenology at the Mz scale. Relaxation of the universality assumptions, 
for instance the one concerning the scalar naassesll3llil, may substantially mod- 
ify the main properties of relic neutralinos ufi-3. We do not discuss these more 
refined supersymmetric models in this paper; only at the end of the next section 
we will compare some results obtained in the MSSM with the ones obtainable 
in theoretical_achemes with and without universality conditions on the soft 
scalar masses t3. 



3 Relic abundance and detection rates in MSSM 

Now we present our main results. Fig.^ shows the neutralino relic abundance 
f2 x /i 2 as a function of m x in the form of a scatter plot obtained by varying 
the parameters over the ranges of Eq. (]?]); £l x h 2 has been calculated with the 
method illustrated in Rem. 

The configurations giving an fl x h 2 in conflict with the cosmological bound 
have already been removed. As expected, large values of the relic abundance 
are provided by configurations with small tan/3 and large rriA, since under 
these circumstances the neutralino pair-annihilation cross section is small. For 
the same reason large tan/3 and small entail small values of fl x li 2 . By 
way of example, in Fig.^j two contour lines denote the values of Vt x h 2 for two 
subsets of neutralino configurations for extreme values of tan/3 and raj. 

We turn now to the presentation of detection rates for two classes of 
searches for relic neutralinos: 
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Direct detection, i.e. measurement of the energy released by a neutralino 
in its scattering off a nucleus in an appropriate detector, by using very 
different experimental technkm^jSLSome of the most recent experimen- 
tal results are given in RefsOa i EjtE3 l E|f3 The rates for direct detection 
are calculated here as explained in Re: 



Detection of the possible signals consisting of up-going muons through a 
neutrino telescope generated by neutrinos produced by pair annihilations 
of neutralinos captured and accumulated inside the Earth and the Sun. 
The evaluation of the muon fluxes, which is a rather elaborate multistep 
process, has been performed here according to the procedure described in 
Ref. El to which we refer for details of the calculations and for references. 



For a presentation of the results concerning a possible excess in the p/p 
ratio in the cosmic rays we refer to Reftj. 

We remind that the detection rates for both direct detection and detection 
of the up-going muon fluxes are proportional to the neutralino local (solar 
neighborhood) density p x . Tcwassign a value to this quantity we employ here 
the following rescaling recipe E3. For each point of the parameter space, we 
take into account the relevant value of the cosmological neutralino relic density. 
When fl x h 2 is larger than a minimal (flh 2 ) m i n , compatible with observational 
data and with large-scale structure calculations, we simply put p x = pi. When 
£l x h 2 turns out to be less than (ilh 2 ) m i n , and then the neutralino may only 
provide a fractional contribution fl x h 2 / (flh 2 ) m i n = £ to Vlh 2 , we take p x = pi£,. 
The value to be assigned to {Q,h 2 ) m i n is somewhat arbitrary, in the range 0.03 < 
(£lh 2 ) min < 0.3. In the present paper we have used (Qh 2 ) min — 0.03. As far 
as the value of pi is concerned, we have taken the representative value pi = 
0.5 GeV • cm" 3 . This corresponds to the central value of a recent determination 
of pi, based on a flattened dark matter distribution and microlensing datfJ: 
Pl = 0.51±g$ GeV • cm" 3 . 

Fig.^ shows a comparison between one of the most stringent experimental 
upper bounds: R e Jf(6 KeV < E ee < 8 KeV) = 1.2 counts/ (Kg day) (this is 
extracted from the experimental data of Ref. til) and the predicted values for 
the same quantity. 

These results are also reported in FigJ^ in a plot versus fl x h 2 , to show the 
expected correlation between the signal and the relic abundance. 



7 



Table 1 : Values of the astro-physical and cosmological parameters relevant to direct detection 
rates. Vr.m.s. denotes the root mean square velocity of the neutralino Maxwellian velocity 
distribution in the halo; V eS c is the neutralino escape velocity and Vq is the velocity of the 
Sun around the galactic centre; p ioc denotes the local dark matter density and (Sl/i 2 ) m i n the 
minimal value of f2h 2 . The values of set I are the median values of the various parameters, the 
values of set II are the extreme values of the parameters which, within the physical ranges, 
provide the lowest estimates of the direct rates (once the supersymmctric parameters are 
fixed) . 





Set I 


Set II 


V r . m . s (km ■ s -1 ) 


270 


245 


V esc (km • s _1 ) 


650 


450 




232 


212 


pioc(GeV • cm -3 ) 


0.5 


0.2 


(17/1 )min 


0.03 


0.3 



In Fig.g we notice that for some configurations, with large tan/3 and small 
rriA, the predicted rates are above the experimental limit. However, we empha- 
size that this feature cannot be used to exclude these neutralino configurations, 
since one has to take into account the large uncertainties affecting some of the 
astrophysical (and cosmological) quantities which enter in the evaluation of the 
detection rates. Table 1 reports two sets of values for the relevant quantities 
of this kind. Set I corresponds essentially to the central values for the vari- 
ous parameters, whereas set II corresponds to those values of the parameters, 
which, within the relevant allowed ranges, provide the lowest detection rates 
(once the supersymmetric variables are fixed). 

For all the results of the present review, unless specified differently, the 
set I has been employed. This is, for instance, the case for the detection rates 
reported in Fig.[| Now it turns out that, if the calculation of the rate -R^f* 
(6 KeV < E ee < 8 KeV) is performed using set II, all predicted rates are 
below the experimental upper bound. However, it has to be stressed that 
the maximum value of the rate falls short of the limit only by a factor of 
2. This obviously implies that even a moderate improvement in experimental 
sensitivities in direct searches may provide an essential information about the 
neutralino parameter space. 

A more direct illustration of this point is provided by an inspection of 
the neutralino-nucleon cross section. Fig.|^ gives a comparison of the predicted 
values for £ • cr^ alar with the relevant upper bound (as derived from the current 
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experimental data Refs0@@@@). 



"scalar stands for the neutralino-nucleon scalar cross-section, responsible 
for coherent effects in processes involving nuclei. Here set I has been used. 
The results of the same calculation with set II for the astrophysical and cos- 
mological parameters are shown in FigJ^. We notice that no predicted values 
are above the experimental upper bound; however, some neutralino configu- 
rations are close to it. We do not consider here explicit upper bounds on the 
spin-dependent neutralino-nucleon cross section as obtainable from present ex- 
perimental data. This point will be discussed in detail in a forthcoming paper 



Let us turn now to the indirect measurements, based on the outcome 
of high-energy neutrinos from neutralino pair annihilation in celestial bodies. 
Here we limit ourselves to the emission from the center of the Earth. The 
flux of the ensuing up-going muons, integrated over muon energies above 1 
GeV and over a cone of half aperture of 30° centered at the nadir, is shown 



m Fig.^ 

10" 14 cm- 



The horizontal 

- x (9o%ai.) 



ie denotes the Baksan upper limit: $ 



Earth 



< 2.1 x 



Also in this case, many configurations provide fluxes in excess of the ex- 
perimental limit. However, if set II of Table 1 is employed, then the number 
of neutralino configurations staying above the experimental bound turns out 
to be marginal. It then follows that present data from neutrino telescopes are 
not yet capable of excluding physical regions of the supersymmetric parameter 
space (contrary to the conclusions of RefEJ). 

Fig.|§| displays $^ ar " 1 versus Rc e to show the expected correlation between 
these two quantities. Also displayed are the locations of configurations which 
would be selected by a model with universality assumptions at Mqut and by 
one where the unification condition on soft scalar masses has been somewhat 
relaxed lil. This illustrates how strong are the constraints imposed by strict 
universality conditions at Mqut and how sensitive are the signals to these 
assumptions. 



From what we have discussed in this part, we can conclude that the present 
sensitivities for the experimental searches examined above are not yet in the 
position of constraining the supersymmetric model. However, the foreseable 
improvements in these sensitivities will make soon possible an exploration of 
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some regions of the neutralino parameter space. 



4 A Higgsino-like relic neutralino 

It has recently been conjectured that the CDF event pp — > e + e~77px@ 
is due to a decay chain involving two neutralino states (the lightest and the 
next-to- lightest ones). The lightest neutralino (that we denote as xakm) has 
been further considered in Refm as a candidate for cold dark matter. If the 
neutralino interpretation of RefEU is correct, then the neutralino would be con- 
fined in a rather narrow region of the supersymmetric parameter space, at least 
as far as some parameters (such as M 2 ,/z, tan/3) are concerned. An important 
question is whether or not a relic xakm could be detected either directly or 
indirectly. Here-we only report a few results of a thorough investigation car- 
ried out in Refill, to which we refer for further details (some detection signals 
f° r Xakm have also been investigated in Refsc3l2il for a limited domain of the 
Xak m parameter space and under the hypothesis that xak m provides a large 
contribution to f2). 

In Ref0 it is shown that for most regions of the parameter space the de- 
tectability of a relic xakm would require quite substantial improvements in 
current experimental sensitivities. Some favorable perspectives for investigat- 
ing a region of the xakm parameter space, around the maximal tan/3 value 
allowed by the model of ReftZl (i.e., tan/3 ~ 2.5), are offered by measurements 
of neutrino fluxes from the center of the Earth and of an excess of p/p in cosmic 
rays. To illustrate this point we display in Fig.^p/p versus $^ art/l for xakm 
configurations with 1.5 < tan/3 < 2.5. 

Diamonds and crosses denote the values of the signals when these are 
within two orders of magnitude from the current value of at least one of the 
relevant experimental upper bounds (at 90 % C.L.). Diamonds (crosses) denote 
values obtained by evaluating the neutralino relic abundance with the exact 
expression (with a low- velocity approximation) . Dots denote the values of the 
signals for the other configurations, calculated in the low-velocity approxima- 
tion only. The experimental bounds are displayed by the horizontal and the 
vertical lines: p/p < 7.5 x 10" 5 E3, ^arth < 2 .1 x 10" 14 cm" 2 s" 1 (90%ai.)^- 

Let us finally remark that the conjecture of Ref@ is certainly very in- 
triguing; if correct, it would have important consequences in particle physics 
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and possibly also in cosmology. However, great caution is in order, because of 
the existence of only a single event.af the -type vp — + e + e~77pT and of the 
non-uniqueness in its interpretationEZrEjO'BE^O. 
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Figure 1: Isomass and isocomposition contour lines for the neutralino in the /i — M2 plane 
for tan (3 = 10. The region excluded by LEP is denoted by dots. m x is in units of GeV. 



13 



Figure 2: Neutralino relic abundance as a function of the neutralino mass. The solid (dashed) 
line delimits the region where neutralino configurations with tan j3 = 55 and m,A = 50-60 
GeV (with tan/3 = 1.1 and rrij\ = 500 GeV) are located. 
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gure 3: flg e (as defined in the text) versus m x . The contour lines have the same meaning 
in Fig.§. 



15 



Figure 4: Rc e versus the relic abundance. 
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Figure 5: The neutralino-nucleon scalar cross section multiplied by the scaling factor £ versus 
m x . Set I of Table 1 is used. The contour lines have the same meaning as in FigM 



Figure 6: Same as in Figj5| except that here set II of Table 1 is used. 
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Figure 7: $E arth versus m x . The contour lines have the same meaning as in FigQ 
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Figure 8: &E arth versus R^e- The solid line delimits the location of configurations for the 
MSSM considered in this paper. The diamonds (the dots) denote the values for a model 
with strict universality assumptions at Mqut (with a deviation from unification condition 
on the soft scalar masses). 
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Figure 9: p/p versus ^^arth £ or -^ AKM configurations. The meaning of the symbols is 
explained in the text. 
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